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Abstract 


The aim of the present work is to introduce a method based on the Cheby- 
shev polynomials for numerical solution of a system of Cauchy type singular 
integral equations of the first kind on a finite segment. Moreover, an esti- 
mation error is computed for the approximate solution. Numerical results 
demonstrate the effectiveness of the proposed method. 
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1 Introduction 


Let us consider a system of singular integral equations of the form 


a) aa fe K(t,r)®(r)dr = F(t), -1<t<1,(1) 


A(t) ®(t) + Bit) | 


where 
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K(t,7) =[Ky(t,7], 45=1,2,-...N, 
FQ) =[flt), fat)... fr @]’, 
&(t) = [¢1(0), G2), ---, on (O)]’, 


Alt) = leax(t)), 4.7 =1,2y00gNy 
Bit) =[by®], t,7 =1,2,...,N. 
Here, (ia pi and {f;}§, are given real-valued Hélder functions and 
{0 es are the unknown functions. The matrices A and B are known such 
that S = A+ B and D= A-— B are nonsingular for all t € [—1,1]. In some 
familiar physical problems, the entries of the matrices A and B are constants. 
The singular integral equations play important roles in physics and the- 
oretical mechanics, particularly in the areas of elasticity, aerodynamics, and 
unsteady aerofoil theory. They are highly effective in solving boundary value 
problems occurring in the theory of functions of a complex variable, potential 
theory, the theory of elasticity, and the theory of fluid mechanics. A general 
theory of the system of equations (1) has given in [12]. 
We study the system (1) in the case that A(t) = 0 and B(t) is a constant 
matrix. Therefore, the ith equation of system (1) takes the form 


N 1 hee N 1 
dy | oO ar + f Kuj(t,7)oj(r)dr = fit), —-1<t<1. (2) 
j=l =a j=? -1 


Studies on this singular integral equation can be found in some literatures 
(see [1,3,6,7]). Chakrabarti and Berghe [3] proposed a method for solving (2), 
using polynomial approximation and collocation points have chosen to be the 
zeros of the Chebyshev polynomials of the first kind for all cases. Kashfi and 
Shahmorad [7] constructed another approximate solution of this equation 
by using the Chebyshev polynomials of the first and second kinds. Some 
other methods for solving this equation can be found in [1,6]. A convergence 
analysis of Galerkin and collocation methods for (2) has been given by Miel 
[11]. 

A special type of (2) is the famous Cauchy singular integral equation 


* 9(7) 


_17T-t 


dr=f(t), -l<t<1l, (3) 


which has the following analytical solutions in four special cases based on 
boundedness of the unknown function ¢ at the endpoints of the interval 
[—1, 1]; see [3,9, 13]. 


Case 1. If the function ¢ is unbounded at the endpoints 7 = +1, then 


_ 4% l * V1-P FH) 
7) = a =. t—T ae eee 
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where do is an arbitrary constant. 


Case 2. If the function @ is bounded at the endpoints 7 = +1, then 


o(r) = eae 


dt , Ler, 


tt = 


and a necessary and sufficient condition of existing this solution is 


Case 3. If the function ¢ is bounded at the endpoint 7 = —1 and unbounded 
at the endpoint 7 = 1, then 


1 = 
ay ae : tLO ay. l<r<l. 
1-7 1+tt—T 


Case 4. If the function ¢ is bounded at the endpoint 7 = 1 and unbounded 
at the endpoint 7 = —1, then 


f1— /l+tf 
ae z dt, 1l<7r<l. 
147 ef. 


An application of (3) were given in [2] by reducing a system of dual integral 
equations to Cauchy type singular integral equations, and more methods for 
solving this equation have given in [5,8, 13,15}. 


In the next section, we investigate approximate solutions for system (1) 
in the above four cases. 


2 Approximate solution 


To find approximate solutions for system (1) in the cases 1,2,3,4, for v € 
{1, 2,3, 4}, we set 


Av ( 
Oj(T) & Gy (7) = MOY Be Palt) FS 1 ce lN (4) 
i=0 


and 
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M 
KE) > OP): 471 Baa Ns (5) 
k=0 
where 
Ti(x) =cos(j),  v=1, 
Uj(0) = Re Y= 
Pyle) = ( ae 
, ___ cos (j+4)0) = 
V; (a) ~~ cos($) ’ v= 3, 
sin((j+4)@) _ 
W; (x) ~~ sin(2) ’ — 4, 


are the Chebyshev polynomials of the first to fourth kinds and 


1, y=1, 
1-7, v=2, 
»v(t) = 1+, vy = 3, 
1—-¢t yv=4 


“) 


in which x = cos(@). The roots of Chebyshev polynomials P,,as4i(a) are 
given by 


tn = (6) 
(2n—1) = 
cos She) , v=, 
2nt _ 
cos (3775) y=A, 
where n = 1,2,...,44 +1. These roots are used as the nodes of Gauss— 


Chebyshev quadrature rules. 


Lemma 1. /10] The Chebyshev polynomials satisfy the orthogonality prop- 
erty 
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0, t#j, 
Tm 1=9=0, v=1, 
» Det 
ih Pii(t)Pyj(t)dt = 4 5, 1=j #0, v=, (7) 
ivl-? 
3 t=), y=2, 
Ty b= 9, v=3,4 


Theorem 1. /10] As the Cauchy principle value of a singular integral, we 
have 


PualT) 4 = i 
7 T= (8) 
Laer Wit), v=3, 


V(t), v=a4. 
Now we describe details of finding approximate solution in cases 1-4. 


Case 1. For v = 1, the relations (4)—(5) take the forms 


M 
#i(1) = vil?) = Gaz DOB MP), F=1,2..N, (9) 
1=0 


and 


Ki;(t,7) := "nl i,j =1,2,...,N, (10) 


where Bj, (7 =1,2,...,N, 1=0,1,...,M4) are unknown coefficients and 
the symbol (3>’) denotes that the first term in the summation is halved. The 
functions 


K,;(t,7)Te( - 
Yaga (t =? ff ; 7 x a, 4,9 =1,2,...,N, B= U)1onacg MM, 


V1—72 


can be determined exactly or may be approximated by using the Gauss— 
Chebyshev quadrature rule, that is, 


2 M-+1 
vik) ~ FAW S> Kig(t, 01,0) Th (21,2) 5 
s=1 


where x1,, obtain from (6). 


36 S. Ahdiaghdam and S. Shahmorad 


Substituting from (9)—(10) into (2) and using (7)-(8) for v = 1, gives the 
system 


N M 1 N M 
SS big By Ui-r( 1 aap et Vign(t je = = flt Ts i=1,2,...,N(11) 
j=1 I=1 j=1k=0 


If the given functions f; and 7,;, are square integrable on [—1, 1] with 


respect to the weight function ae then they can be expanded as 


vige(t) & D9 | Gages Ui(2), t, § =1,2,...,N,k=0,1,...,M, 
(12) 
14(t)~ Mp cat), i=1,2,...,N, 


where the coefficients 
Giga = 2 [2 VI— Praga (t)Ui(t)at 


Ae fey fea yf HES Kay (t, TU (t) Ta (7) drt 
1912 Ne. SOU Ag M1 RaSO No 1, 


I 


ty = ee a V1 —# fi(f)U(t)dt, i=1,2,...,N,1=0,1,...,M@—-1, 
can be approximately determined from 


Giget ~ Certara Dosen (1 — 03,,) Kig (@2,r 81,9) U1 (@2,r) Tk (1,5); 
(13) 
M 
Cal © aapgay Dorai (1 — £3.) fi(e2,r)Ui (2,1). 


Using (12) in (11), we have 


M-1N 4 N M ; M-1 
S- S © big By qu4a3Ui(t) + 5 a aes Giz BjrUi(t) > Cil Ui(t 
l=0 j=l 1=0 j=1k=0 


which leads to the linear system 


N M 
1 ’ F 
> bi Byitiy +5 D0 Gijri Pik =cyj,1=1,2,...,N,1=0,1,...,.M—-1, 
j=l k=0 
(14) 


for the unknown values 8j, (7 = 1,2,...,N, k = 0,1,...,M). Tak- 
ing $11,...,8N1 arbitrary, the remaining coefficients 6;, are uniquely found 
from the linear system (14) which determine the elements of the vector func- 
tion ® via (9). 
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Case 2. We set v = 2 in (4)-(5) and substitute them in (2) to get 


N M 1 N M 1 
SOS bi BT) +5 dodo walt ein = — ~filt), i= Ly 2iesegt¥ (15) 
j=1 1l=0 j=1k=0 


where we used the formulas (7)-(8). Then, we expand the functions f; and 
Vijk as 


yun (t) & oy Gum T(t), 1,9=1,2,...,N, k=0,1,...,M 
1 M'! . 
iat) ay CHI 1=1,2,...,N, 

where the coefficients are determined by 


Giga = 2 f°, ere) at, 4,7 =1,2,...,N, kL =0,1,...,M, 
= 4% ft fh fe Kay (t, 7) Us (7) drat 


1-t? 
a= 3 f' seph(OT@ dt, 1=1,2,...,N, 1=0,1...,M, 
or approximated by 


M+1 M-+1 
Gijki & WEDD r= 7 ; Gel — x3 s)Kig(21,r, €2,0)Ti(£1,r) Ug (22,5), 


M+1 
ca sag Dra fi(z1,r)Ti(21,r). 


Using the last expansions in (15), returns the following linear system of equa- 
tions 


1 N M : 
3 Loja1 ko Gazi Bjk = Cit, 42 12.2, 51=0, 


ae {—bi;8;-1 1 § hey GismsBje} = Cil; i= 1,...,N, l= 1,...,M, 


for the unknown values 6;, (j = 1,2,...,N, 1=0,1,...,M). Then ele- 
ments of the vector function ®(¢) obtain from (4). 


Cases 3,4. Similar to cases 1 and 2, we get the linear systems 


N M 
{but +3 Gant =e i=1,2,...,N, 1=0,1,...,M, (16) 


j=l k=0 


and 
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N M 
3 { but +3 Gant =n i=1,2,...,N, [=0,1,...,M,(17) 
j=l k=0 


respectively for v = 3 and v = 4, and then we determine the elements of 
corresponding vector ® via (4). 


3 An estimation error and numerical results 


In this section, we describe an estimation error for the approximate solution. 
Let 


B(t) = [yi l(t), pelt), .--, eu (O)]7 


be the vector of approximate solution of the system (1) and let E(t) = ®(t) — 
®(t) be the associated vector valued error function. Due to the approximation 
®(t), for A(t) = 0, system (1) may be written as 


B(t) i. 20) dr + [. K(t,r)®(r)dr = F(t) + H(t), -1<t<1, (18) 


where the perturbation term H(t) obtains from 


H(t) = Be) f . 


o) dr + [. K(t,7)®(r)dr — F(t), -1<t<1. 


TA 


Subtracting (18) from (1) yields a system of error equations as 


B(t) 1s 0) dr + [. K(t,r)E(r)dr = H(t), —1<t<1, 


TA 


which is solvable approximately like system (1). 


The following examples illustrate application of the method. 


Example 1. Let 


A(t) = 0, Bet) =| qa]. xen= |" él, Feo= |"). 


and find the solution of system (1) in case 1. 


By the above information, the system (1) is reduced to 


ia — dr + fic —t)d1(r)dr + fo, th2(r)dr =7, —-l<t<l, 
(19) 


fy ea dr + ‘i Th1(T)dt + fic + t)ho(r)dr = Int, -1<t <1, 
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and since the matrices S = A+B = Iz and D = A—B = —J, are nonsingular, 
then system (19) has a unique solution. The kernels K1,(t,7), K2;(t,7) 
(j = 1,2), and the functions f; and f2 are polynomials of degree at most 1, 
so we set 


ds(1) = Fs {BjoTH(7) + BTU) +8 T)}, §=12, (0) 


and 
Kij(t,7) = yaj0(t)To(7) + i t)Ti(7), i, 9 = 1,2, 
F(t) = qoVo(t) +. cnt), i= 1,2, 
where 
yu1o(t) = —ZUi(t), Yx11(t) = Uo(t), Y20(t) = quilt), yi2i(t) = 0, 
y210(t) = 0, yoi1(t) = Uo(t), ya20(t) = 5Ui(t), yea1( 


Substituting these expansions into (19) and using (7)-(8), for v = 1, we 
obtain 


By Uo(t) + Br2U1(t) — 5 810U1 (t) + §811Uo(t) + §520U1 (t) = Volt), 


B21 Up(t) + B22Ui(t) + $811U0(t) + $820U1(t) + $B21U0(t) = Ui(t). 
Then the linear independency of {Uo(t), Ui(t)} implies 


3814 — 1, 

—4B10 + Bi2 + $820 = 0, 
3511 + $821 = 0, 

5520 + B22 = 1. 


A nonunique solution of this system for the arbitrary values of 319 and Sg 
is given by 


= Biz = 5 (B10 ~ B20) B20, Pie ey ee 


Bio, Buz 3 F 5 


For example, if B19 = G20 = 2, then 812 = 622 = 0 and so we find from (20) 
that 


2742 —274+2 
oi(T) = > ba(T) = =: 


(see Figure 1 for the behavior of these solutions). 


Example 2. Solve the problem of Example 1 in the case 3. 
In this case, we set. 
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Figure 1: The plots of approximate solutions of Example 1 for M=2. 


by(7) =f 52 (Bo Volr) + BE, F= 12, 21) 
and 

Kij(t,7) = vigo(t)Vo(7) + Vi QVi(7), 7,9 = 1,2, 

fi(t) _ cio Wo(t) + ciiWi(t), i= 1,2, 
where 
‘ya1o(t) = Wo(t) — sWilt), — -Ya11(t) = Ye10(t) = Yo11(t) = Ya21(t) = 5Wolt), 
yi20(t) = —3Wo(t) + $Wi(t), Y21(t) =0, —Yaa0(t) = 5 Wi(t), 
cio(t) = 1, ci1(t) = 0, C20(t) = —l, co1(t) =i. 


Substituting these expansions into (19) and using (7)—(8) for v = 3, result 


BioWo(t) + B11Wi(t) + BioWo(t) — 5 F10Wi (t) 
+5811Wo(t) — 5820Wo(t) + 5820Wi(t) = Wo(t), 


B29Wo(t) + B21Wi(t) + $810Wo(t) + $611Wo(t) 
+5 B20W1 (t) + 5821Wo(t) = —Wo(t) + Wilt), 


and from the linear independency of {Wo(t) and W;(t)}, we get the algebraic 
system 


2810 + 5811 2 320 =1, 

~ 3510 + B+ 720 = 0, 
$10 + $811 + B20 + $621 = -1, 
3820 + B21 = 1 


oI 


which has the unique solution 
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Figure 2: The plots of approximate solutions of Example 2 for M=2. 


10 28 22 20 


Pio = — 5 Bu = 57> Po =—> Pu => 


and the solutions of (19) can be found via (21). The graphs of these solutions 
plotted in Figure 2. 


Example 3. Solve the following system in the case 3 [16]: 


a oul7) dr+f- da(7) dr = (t?7 +1)sin2t, -1<t<1l, 


see -1 T- 
eae Ou) dr 4 ie ga(7) dr = tcos 2t, -1<t<1l. 


In this case, we will approximate the unknown functions as 
M 

bi(t) = pit) = VI 7?) Big Uj(7), 1 = 1,2, 
j=0 


where the exact solutions are 


—7r2 rl 241) sin2 cos 2 
dur) = — VG pt, Ca sinatteos2t ap, 


= vice 3t cos 2t—2(t 241) sin 2t 
2(T) = ie Foe > dt. 


We define the error functions as 
E,(r) = |¢i(r) — gi(r)|, 2 = 1,2, 


where the exact solutions ¢; and @2 are calculated by using the Maple code 
int (expression,x=a..b, CauchyPrincipalValue). 
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Table 1: Comparison of our results (E;) with the results of [16] (¢;) for M=8. 


T p1(T) $2(T) Ex(r)  Ea(r) e1(r) — e2(r) 
—0.96 —0.2164694659571311 0.5107921932649533 le—8 2e—8 8e—5 2e—4 
—0.70 —0.5394637532569845 1.185306580331174 3e—7 6e—7 4de—5 le—4 
—0.25 —0.5831756914063946 1.129312022840245 4de—7 9e—7 9e—6 2e—5 

0.25 —0.5831756914063946 1.129312022840245 4de—7 9e—7 9e—6 2e—5 
0:75 —0.5394637532569845 1.185306580331174 3e—7 6e—7 4de—5 le—4 
0.96 —0.2164694659571311 0.5107921932649533 le—8 2e—8 8e—5 2e—4 


T pi(T) p2(T) Ey (rT) Eo(r)  e1(T) — €2(T) 


—0.96 —0.2164694659571311 0.5107921932649527 0 6e—16 - - 
—0.70 —0.5394637532569844 1.185306580331173 le—16 le—15 - - 
—0.25 —0.5831756914063935 1.129312022840244 le—15 le—15 - - 


0.25 —0.5831756914063935 1.129312022840244 le—15 le—15 - - 
0.75 —0.5394637532569844 1.185306580331173 le—16 le—15 - - 
0.96 —0.2164694659571311  0.5107921932649527 0 6e—16 - - 


In Table 1, we compared our numerical results (absolute errors E, and 
E2) with those reported in [16] (absolute errors ¢; and €2). Table 2 shows 
our results for M = 15, which were not reported in [16]. 


Example 4. Consider the singular integral equation 


1 
‘ | or 
dr +f e'd(r) dr = 1-27 +i + xe -l<t<l, (22) 
1 


= rat 


with the exact solution ¢(7) = ¥ ior (27 — i) in the complex plane, where 
i= /-1. 

By taking ¢; := Re{¢} and ¢2 := Im{¢}, equation (22) is reduced to the 
system of singular integral equations 


f2, BO art f*, cos(t)o1(7) dr — f7, sin(t)da(r) dr = 1-242 + 4sint, -1<t<1, 


fn eek) dr + fet, sin(t)¢1(7) dr + ia cos(t)@2(7) dr = t — 4 cost, -l<t<1. 


Similar to Examples 1 and 2, by setting 
M 
bi(t) = pit) = VI 7?) [Big Uj(7), 1 = 1,2, 
j=0 


we get the exact solutions for M = 1. 


Example 5. Consider the system of singular integral equations [14] 
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1000 1, BO dr + 10 fl, SO dr = fit) +in(t), -1<t<1, 


(23) 


500 ft, 2 dr + 200 ph 2) dr = fo(t) + igo(t), 1 <t <1, 


where 


53511 t 
4000 


53929 


fi(t) = —990¢8 + 1089¢7 + 937¢° 2000 ° 


26704 4.5 34916144 | 79232743 1761 42 
25 t 1000 et 2000 ot 250 t 


30873 42 695334 4 
Toooo & + Zo00 ! ' 


1130501 


7 
1189¢ 40000 ° 


_ 8 897146 , 11904745 , 16396144 _ 279198 ,3 
gi(t) = 990¢ To & + “too! + “Saat 6a5 


= 8 7 6 2447,5 _ 8607;4 , 73543 _ 17253,2 4 29541 14701 
f(t) 300t° + 330¢° + 215¢ tot — oot + yet 3000 & + “4000 & — “to00> 

- 8 7 6 | 1462,5 , 9419;4 _ 2754943 _ 18342 7957, 1 57583 
g2(t) = 300t° — 380t" — 197” + Pet? + Aap t 350 — aool — ao00! + “ao00 - 


It is easy to see that system (23) is equivalent to the following disjointed 
singular integral equations, 


Lf MOEN a = ROSH, 1 cre 
oe Im{¢1(7)} [a 20gu() ga) teppei 
nT JI Tt ? , 

24 
L pl, Refda) gp = 2-H) _y ct <1, (24) 
2) ee ge Poa ete, 


The exact solution of system (23) was reported in [14] for the case r = 4. 
By applying our method for the equation in (24), we get the real parts of 
unknown functions exactly but, for the imaginary parts, we obtain the error 
functions 


1l-—T . 
E;(r) = |Im{¢gi(r)} — Im{yi(r)}| = rE hi(r), 1=1,2, 
a ad 
where 
hi(r) = 5.128205128205128 x 10-8(1 +7), -l<r<l, 
ho(T) = 5.128205128205128 x 10-®(1 +7), -l<r<l. 


The plots of hy and hz and the regular parts of the error functions, are 
shown in Figure 3. As shown in these plots, the absolute errors increase 
when 7 approaches to 1, because in case 4, the singularity for the integral 
part happens at 7 = 1. 


Example 6. Consider the problem of a half plane containing a crack parallel 
to the boundary, which is illustrated in Figure 4 and formulated as the system 
[4] 
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The plot of h_1(t) The plot of h_2(t) 
1.x 10°74 0.000010 
8.x 10°84 0.000008 
6.x 10784 0.000006 
4.x 10°54 0.000044 
2. x 10°84 0.0000024 
0 : 0 
1 0.5 0 0.5 1 -1 05 0 0.5 1 
7 t 


Figure 3: The plots of the regular parts of the error functions in Example 5. 


+, 2 drt 2, [Kis(t,7)b1(7) + Kia(t,7)b2(r)dr] = 0, 
(25) 


* 22 art fi, | Kart, r)¢1(7) + Koalt, 7) 62(7)dr] = 7, 


with 


_ 8h? (rt) 4h? (r—t)[12h? —(r—t)?| 
Kir (t,7) a (rT DEE? + [(r— D2 -4h2 2 (Caza ; 


8h3 [4h? —3(7-t)’] 


Kyolt, T) = Koi (t,7) = [(7-t)24+4h2)>? 


) _ enue 8h? (r—t) 4h? (r—t)[12h?—(r—-t)?] 


Koa(t,7) = — Gaya — [GaP (oy oan 


where fh is the distance of crack from the boundary. The physical conditions 
of the problem impose that the relations 


i: gi(t) dr =0, i g2(T) dr = 0, (26) 


gilt) = bi(—t), G2(t) = —ba(—t) 


are satisfied. Therefore the unknown functions may be expressed as 


and 


1 1 
oi(T) = Fee Sata » a(t) ~ a i(7). (27) 


FT] 
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Figure 4: Crack parallel to a free boundary in Example 6. 


For v = 1, it follows from the orthogonality condition (7) that the second 
condition in (26) satisfies and the first one gives 319 = 0. 


Taking 611 = 821 = 0, the remaining coefficients 6;; are uniquely deter- 
mined from the linear algebraic system (14) for each values of h and M. This 
leads to find the functions ¢; and ¢2 from (27). 


The stress intensity factors 


ky = lim V1-— 7? do(r) 


T1- 


kg = lim 1-7? ¢,(rT) 
tT O17 


and their absolute estimation errors(Est.Err.) are reported in Table 3. For 
h = oo and K;,;(t,7) = 0, from (27) and (13)—(14), the exact solutions of (25) 


are obtained as - 


T) =0, T) = ———_, 
or(7) 2(T) a= 
which give kj = 1 and ky = 0. This is shown in the last row of Table 3. 
The table shows the rapid convergence of the results even for relatively small 
values of M. 


4 Conclusions 


We described a new idea of using Chebyshev polynomials for the numerical 
solution of the system of singular integral equations of the first kind. In 
Section 3, we illustrated this idea by using system of different kind of singular 
integral equations (Examples 1-5). In Example 6, we studied a crack problem 
in solid mechanics and reported the numerical results (see Table 3) to show 
the efficiency and rapid convergence of the proposed method for all these 
kinds of problems. 
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Table 3: Stress intensity factors for the crack parallel to the boundary 


h M ky Est.Err. ky, ko Est.Err. ko 
0.2 6 4.878800637605022 6.3e-14 1.750099102171126 6.2e-14 
7 4.788277537335018 1.1e-14 1.727809740547429 4.1e-15 
8 4.760729834685963 4.8e-14 1.719782910590219 1.1e-14 
0.4 3 2.607272141646415 4.3e-15 0.7745787927510580 1.6e-16 
4 2.594500911475041 7.3e-15 0.7266641783709941 5.0e-15 
6 2.594423234973139 4.2¢e-14 0.7376171346942053 3.3e-16 
0.6 2 1.834057544899021 1.1e-15 0.5664257041432605 4.1e-16 
5 1.960455689663461 6.1e-15 0.4297949760368867 1.6e-15 
0.8 2 1.608371955353828 1.6e-15 0.3323260582700188 2.8e-16 
3 1.660617572058080 8.7e-16 0.2675691556476836 4.6e-16 
1.0 2 1.461157081431933 2.0e-15 0.2104682299562445 1.1e-16 
4 1.485914720666516 2.1e-16 0.1796691052492212 1.1e-16 
1.2 4 1.372176156193755 5.0e-16 0.1234414146531335 0. 
1.5 4 1.262800608570183 1.6¢e-15 0.07465158121522054 1.7e-16 
2.0 3 1.162112249974693 1.le-15 0.03662808437088003 0. 
3.0 2 1.077621553329114 3.3e-16 0.01274529646673066 0. 
10 2 1.007451045420713 2.6e-16 0.00037197952964307 0. 
oo i 1 0 0 0 
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